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Abstract. Using a tight-binding atomistic simulation, we simulate the recent atomic-force microscopy experiments 
probing the slipperiness of graphene flakes made slide against a graphite surface. Compared to previous theoretical 
models, where the flake was assumed to be geometrically perfect and rigid, while the substrate is represented by a static 
periodic potential, our fully-atomistic model includes quantum mechanics with the chemistry of bond breaking and 
bond formation, and the flexibility of the flake. These realistic features, include in particular the crucial role of the flake 
rotation in determining the static friction, in qualitative agreement with experimental observations. 

PACS. 6 8.35.Af, 62.20.Qp, 81.05.Uw, 07.79.Sp 



1 Introduction 

The scanning tunneling microscope (STM) |[T], and even more 
the atomic -force microscope (AFM) ||2|, have triggered per- 
haps the biggest wave of advances and discoveries ever in sur- 
face science and nanoscience. Experimental investigations of 
friction on the atomic scale have become possible by virtue of 
the friction force microscope (FFM). In a FFM experiment a 
sharp tip scans a sample surface with atomic precision, while 
lateral forces are recorded with a resolution that can reach the 
pN range. Since Mate et al. [3] investigated the nanoscale peri- 
odic frictional force map of a graphite surface using a tungsten 
tip, many studies have been conduced experimentally and the- 
oretically. In recent works the Leiden group II4I5I has probed 
quantitatively the well known slipperiness of graphite, respon- 
sible for its excellent lubrication properties. Morita et al. ||6l 
suggested that in FFM experiments on layered materials, such 
as M0S2 or graphite, a flake, consisting of several hundred 
atoms in contact with the substrate, can attach to the tip. By 
controlling the relative angles of individual nanoflakes to achieve 
a suitable lattice mismatch, thus incommensurate contact [7], 
an almost frictionless sliding was demonstrated for dry and 
wearless tip-surface contact, a phenomenon known as super- 
lubricity. Several experiments and calculations an have been 



ulated FFM experiment where dissipation of a finite graphene 
flake is driven along a graphite substrate, in a fully atomistic 
scheme based on a tight-binding (TB) force field. Compared to 
similar models in the literature 1161171181191 . where the flake 
is assumed to be perfect and rigid and the substrate is rep- 
resented by an analytically defined static periodic potential, 
our fully-atomistic TB simulation explores two realistic fea- 
tures, namely: (i) The flake-substrate interaction potential is not 
classical and contains quantum mechanics with the chemistry 
of bond breaking and bond formation, (ii) The flake is non- 
rigid, so that during its advancement it can deform and relax. 
In Sec.|2]we introduce the model implementation details. Sec- 
tion [3] reports the results obtained, in particular for the friction 
dependency on the flake size, the rotation angle relative to the 
substrate, and the applied load; we compare the results of the 
present model to experiment and to previous calculations. The 
final Sec.|4]discusses the results and the advantages and draw- 
backs of the present model. 



2 The model 

We describe the sliding by means of a generalized Tomlinson- 
like model similar to that of Ref. ITTll . but including the fol- 



probing the effects of lattice mismatch on friction 181911011 II 1211 Sl owing features; (i) the interaction among all carbon atoms is 



showing that incommensuracy often prevents global ingrain- 
ing of large areas, thus attenuating the consequent strongly dis- 
sipative stick-slip motion. Theoretically, atomic-scale friction 
on ideal solid surfaces is often described by simple balls-and- 
springs models such as the Tomlinson model lll4l . where a sin- 
gle atom, or a more structured tip ifTsl . is dragged through a 
spring along a static periodic potential energy surface. In the 
present work we implement a Tomlinson-Iike model of a sim- 



realized in terms of the tight-binding scheme of Xu et al. Il20l . 
and (ii) the flake is not rigid but is allowed to deform and 
rotate while sliding. Interatomic forces are computed as cus- 
tomary in the TB scheme IISTI : the hopping parameters and 
the pairwise repulsive potential term follow the scaling form 
given by Xu et al. 1201 . AH interatomic interactions vanish at 
a cutoff distance = 2.6 A. This distance sits in between 
the nearest-neighbor and the next-nearest-neighbor distances of 
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Fig. 1. (Color online) The definition of the angle (j) measuring the ini- 
tial rotation of the flake (thick, pink lines) with respect to the substrate 
honeycomb structure (thin, black lines). The pulling line is defined by 
its distance h' to the center of a substrate hexagon and the pulling an- 
gle 6. To take as a reference the AB stacking (0 = 0), we define h' in 
terms of the parameter h = h' + \dgrij,i,. 
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Fig. 2. (Color online) The computed normal force per atom 

atom 'IS a 

function of the fixed rigid flake-substrate distance d. The curves refer 
to different stackings (AA or eclipsed, AB and an incommensurate 
obtained starting from AB rotated by (p = 30°) for a 96-atom flake 
undeformed flake. The compression curve based on the low-pressure 
data of Ref. II26I is reported for comparison. 



carbon atoms of the equilibrium sp-' diamond structure, and of 
the sp- graphene plane. It is also shorter than the interlayer dis- 
tance of graphite, which is as long as 3.35 A 1221 . The present 
TB model has been applied successfully to investigate several 
low-dimensional carbon systems 1231241251 . In particular, this 
parameterization reproduces the experimental equilibrium dis- 
tance dgraph = 1.4224 A of the graphitic plane. To study fric- 
tion, we use a model constituted by a graphene flake sliding 
over a single infinite rigid graphene sheet. The infiniteness of 
the substrate is simulated by repeating periodically a regular ar- 
rangement of A^^"*" atoms at positions r?"'' in the x-y plane. We 
consider a periodic rectangular supercell, containing as many 
carbon atoms as necessary for the flake not to interact with its 
periodic images. 

We have studied three different regular hexagonal flakes 
composed of A^" = 24, 54, and 96 atoms, respectively. Initially, 
the flake is rotated by an angle (p with respect to the substrate 
and translated horizontally to put its center of mass along a slid- 
ing line at a distance li' from the center of a substrate hexagon. 
As illustrated in Fig. [T] this line is oriented at an angle 6 from 
the X direction, which is defined by being parallel to the zig-zag 
direction of the honeycomb lattice. In our simulations we drag 
the flake along this pulling line and usually let the flake atoms 
relax in all directions {x, y and z), whereas the substrate remains 
completely rigid. The angle defines the stacking mismatch, 
which has a central importance for friction. Specific values of 
6 such as 0° and 30° define special pulling directions where the 
pulled flake encounters periodic corrugations, while aperiodic 
corrugations are experienced for generic 6 angles. The range of 
interest of both angles and (p runs from 0° to 30°: outside this 
range we recover equivalent geometries. 



To simulate an AFM experiment we introduce a constant 
load Fff pushing the flake against the substrate along the verti- 
cal z direction and simulating the force applied by the actual tip. 
This load acts against the reaction forces produced by the inter- 
action with the substrate. These forces are reported in Fig.|2]for 
a A^" = 96-atom flake in several configurations, and of course 
vanish beyond = 2.6 A, the TB interaction cutoff. For a 
distance d > 0.18 nm, the load force per atom does not ex- 
ceed 10 nN (a total load in the tens to few hundred nN for a 
flake composed by 10 to 10^ atoms, corresponding to a load 
pressure ^ 4 Mbar), which is a value practically accessible to 
FFM experiments 08191101 . A load per atom of 0.5 nN, withing 
the selected force-field model, produces an approach distance 
near 0.21 nm, similar to the one obtained by assuming for the 
flake-substrate system the equilibrium interlayer separation of 
graphite (3.35 A) and the very soft c-axis compressibility 1261 . 
expressed as dlnc/dP - -2.8 x 10"^ bar"'. Note however that 
we apply this compressibility relation, also sketched in Fig. |2] 
for distances and pressures that go beyond its linear-response 
range of validity: in the region of close approach, d < 0.22 nm, 
the actual force is likely to increase more rapidly, like in the 
TB model curves. We implement a Tomlinson-like dynamics 
with each flake atom coupled horizontally to a rigid "support" 
by elastic springs. The support is a set of ideal graphene-net 
points, which coincide with the initial flake atomic positions: 
j.^^"P _ j.n^f _ Qy -pjjjg support is then translated rigidly parallel 
to the substrate. Its orientation is fixed once and for all by the 
angle (p. The support advances by steps of length Sx along the 
direction defined by the pulling angle and lateral shift h'. Af- 
ter a few tests, we select an advancement step 6x = 0.0024 nm. 
After each advancement step relaxation, we evaluate and store 
the total spring energy and the total dragging force, defined as 



Federico Bonelli et al.: Atomistic simulations of the sliding friction of graphene flakes 



3 



follows: 
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where the // symbol indicates the in-plane component. The 
component F^^^ of F'j^'^ along the pulling direction equals the 
force needed to make the support advance, so that the work of 
this force equals F'^j^^ 6x, for an infinitesimal displacement in 
the pulling direction. 

In an AFM experiment, the scanning tip speed is typically 
of the order of tens or hundreds nm/s, much slower than the fast 
dynamics of the flake. Assuming that the substrate temperature 
is fairly low, it is appropriate to consider a quasi-static motion 
of the flake as follows: after each advancement of the support, 
all flake atomic positions are made relax in all directions by 
damped dynamics Q under the combined action of (i) the TB 
forces, (ii) the vertical load force Ff^, and (iii) the horizontal 
spring forces that attract the flake atoms near the support points. 
When the stationary equilibrium position is reached (defined by 
no force component exceeding a threshold of 10 ~ nN), the sup- 
port moves one step further and the whole relaxation procedure 
is repeated until the support reaches the end of a pre-defined 
path. With the selected fairly small advancement step 6x, each 
relaxation requires typically 10 to 200 MD steps. We consider 
a path of moderate length d ^ I nm divided into approximately 
400 advancement steps. The spring elastic constant coupling 
the support and the flake in the x-y plane is a quite critical pa- 
rameter of the present model. Softer springs allow the flake 
a greater freedom to translate, rotate, and deform, with better 
pinning to energetically favorable sites and more pronounced 
stick-slip motion and higher friction. Harder springs enforce a 
more stiff flake showing little or no stick-slip motion. The limit 
— > oo matches the model by Verhoeven et al. ifTTl . The spring 
constant mimic the combined interaction between the flake and 
the tip and the elastic tip response. We suggest that a value 



eV N 
8.01 - 

A2 m 



(4) 



corresponding to about 10% of the stretching stiff'ness of a 
carbon-carbon bond within a graphene layer should probably 
be fairly realistic. We also perform computations with softer 
{K = 0.\ eV/A^) and harder {K, = K, = 2.5 eV/A^) springs, as 
detailed in Sec. [3] 

In our calculations, by convention, we estimate the static 
friction force Ff^c along a given sliding path by the force level 
below which 95% of the spring force values Fj^ encountered 



' The advantage of damped dynamics with respect to a standard 
energy-minimization algorithm is that at each step it relaxes smoothly 
and predictably to the nearest local-minimum configuration. On the 
other hand, this algorithm is computationally less efficient than, e.g., 
a conjugated-gradient technique. Eventually we settle on a fairly fast 
converging damped dynamics with a time step St = 4 fs. 




tip position x^.^ [nm] 

Fig. 3. (Color online) The spring total lateral force component F^^""^ 
projected along the scan line, for a 24-atom flake dragged starting 
from an initial stacking AB along a pulling angle 9 = 15°, and with 
a stacking angle (p = 30°. The spring constant is = 0.1 eV/A^, the 
total load Fn = 100 nN. The 95% force level (dashed line) estimates 
the static friction force Ff,ic. 



forward scan 
backward scan 




tip position x^.^ [nm] 

Fig. 4. (Color online) The dynamic-friction force is evaluated as the 
energy dissipated through a forward-backward scan (solid and dashed 
lines respectively). The shaded area between the two curves measures 
the energy dissipated by friction. The conditions are as follows: A'" = 
24, Fn = 100 nN, 6* = 0°, = 0°, K = OA eV/A^ 



along the path (at the end of each relaxation, and excluding 
an initial transient) are found. This definition is illustrated in 
Fig. |3] The dynamic-friction force is defined in terms of the 
energy dissipated in a forward-backward loop. Figure |4] illus- 
trates this concept: the component of the spring force parallel 
to the advancement direction shows a clear hysteretic behav- 
ior in a forward-backward scan. The shaded area enclosed be- 
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tween the two curves measures the energy Efn^ dissipated by 
friction, and is clearly related to the stick-slip events. The av- 
erage dynamical friction force Fj'J^"''™'^ = Efnc/d ^ 0.325 nN 
is of course smaller than the static friction force Ffrfc - 0.9 nN 
evaluated according to the 95% protocol described above. In 
the following we will focus on the static friction force Ffric, 
which is cheaper to compute and eventually of the same order 
as F'. 



dynamic 



. As the flake advances along its path, it deforms and 
rotates around its center of mass. In particular, to understand 
the evolution of the static friction force, it is useful to track the 
flake instantaneous stacking angle 0^, which generally differs 
from the fixed support angle 0. At each relaxed configuration, 
we calculate (pA as an average over all flake atoms ;, of a sine 
projection obtained as the length of the vector product of two 
vectors in the horizontal xy plane: R™, joining the flake cen- 
ter of mass to the ;-th flake atom, and the corresponding vector 
R^}™ computed for the unrotated = 0° support. Explicitly: 



arc sin 



pOcm X 



_ ^cm xg 



,0 cm V \ 



(5) 



3 Results 



We analyze the friction force dependence on several physical 
parameters, namely the pulling angle 9, stacking angle 0, ap- 
plied load, flake size and position of the scan-line. Experimen- 
tally, the number of flake atoms, estimated in the order of 100, 
is not well determined, while the total applied load and the total 
force acting on the flake are under control in the FFM. For ease 
of comparison with experiments, our discussion shall always 
deal with total quantities, i.e. summed over the flake atoms. 



3.1 Relaxation to equilibrium 

Optimally stacked configurations are important in providing 
the most efficient sticking points during a friction sliding ex- 
periment. Figs.|5]|6] and|2]report typical such relaxed configu- 
rations for 24-, 54- and 96-atom flakes respectively, obtained 
under the action of a total load of 100 nN, and with no spring 
connection to a support {K - 0). The relaxed configuration 
(up to a symmetry rotation/translation) depends only moder- 
ately on the starting stacking, unless the initial stacking angle is 
strongly incommensurate. The average vertical flake-substrate 
separation is 0. 192 nm. The equilibrium configuration tends to 
arrange the flake so as to minimize the number of flake atoms 
stacked on top of a substrate atom: indeed Fig.|2]shows that the 
ecHpsed "AA" stacking produces the strongest repulsive force 
at the same distance. Geometrically different non-optimal con- 
figurations are characterized by typical excess total energies of 
1 eV or less. 



3.2 The stick-slip movement 

We come now to the actual simulation of sliding friction: to 
start with, Fig.[8]displays the pulhng force Fj' measured along 
two sliding paths of different commensurability nature: = 0° 




Fig. 5. (Color online) The relaxed configuration of the 24-atom flake. 
Static substrate atoms are clear (white), flake atoms are dark (red). 




Fig. 6. (Color online) The relaxed configuration of the 54-atom flake. 
Static substrate atoms are clear (white), flake atoms are dark (red). 



(periodic) and 6 - 15° (aperiodic). As expected, the regular 
stick-slip pattern of the 9 = 0° path is replaced by an irregular 
dependency in the 9 = 15° trajectory. The initial part of the 
9 = 0° trajectory is not periodic because of the usual startup 
transient behavior: the first 0.2 nm are omitted from the calcu- 
lation of the friction force, as discussed above. 

The stick-slip motion and the role of the spring stiffiiess 
is understood even better by comparing other physical quanti- 
ties with FJ\ In particular, Fig.|9]displays the internal energy 
of the flake-substrate interaction, the displacement of the flake 
center of mass along the pulling direction, and the actual stack- 
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Fig. 7. (Color online) The relaxed configuration of the 96-atom flake. 
Static substrate atoms are clear (white), flake atoms are dark (red). 




tip position x ^.^ [nm] 

Fig. 8. (Color online) The component of the springs force Fj"^ in the 
dragging direction as a function of the support advancement distance 
Xtip for two different pulling directions: 6 = 0° (solid) and 6 = 15° 
(dashed). The simulation involves a 24 atom flake with total applied 
load of 100 nN, support stacking angle = 0° and spring constant 
= 0.1 eV/A^. 



ing angle (pA- We focus initially on the solid curves, obtained in 
a simulation based on soft springs with K = 0.1 eV/A^. After 
the initial transient, where the flake explores once a configura- 
tion with a negative <pA, it then jumps back and forth between 
two kinds of sticking configurations: the most favorable one 
characterized by 0^ ^ 4°, and another one, at 0.1 to 0.2 eV 
higher energy, near (pA - 0°. These stick-slip jumps overcome 
energy barriers whose heights are of order 0.2 to 0.3 eV. This 
energy amplitude sets the temperature range of validity of the 
present zero-temperature calculations to a few hundred Kelvin: 
when thermally-activated slips through energy barriers do not 
have enough time to occur, i.e. for a not-too-small tip advance- 
ment speed 1271281 . these slips are unlikely and our estimates 




0.5 
tip position x 



up 



Fig. 9. (Color online) Several physical quantities plotted as functions 
of the support position j,ip: (a) the parallel component of the springs 
lateral force Fj\ (b) the flake excess energy - is^f' - E'^'*, (c) 
the flake center-of-mass advancement xcm along the pulling direction, 
and (d) the flake instantaneous rotation angle (p^ relative to the sub- 
strate. Two different values of the spring constant are compared: soft, 
a: = 0.1 eV/A^ (solid, squares), and hard, K = 2.5 eV/A' (dashed, 
triangles). These simulations involve A'" = 24, = 100 nN, 9 = 0°, 
and = 0°. 



of friction should be fairly reliable. Stick-slip events occur with 
correlated jumps in the spring force, flake excess energy, flake 
position, and stacking angle. The very different dashed curves 
show that stiff springs (with K - 2.5 eV/A^) produce a much 
stronger and more rigid binding of the flake to the rigid sup- 
port. Accordingly, such an unrealistically rigid coupling sup- 
presses the stick-slip behavior: both the advancement, shown 
in Fig. HJc), and the spring force. Fig. |9ja), become smooth 
and jump-less. Despite the suppression of stick-slip, we ob- 
serve higher force peaks for the stiffer springs, thus indicating 
a higher static friction than for the softer springs. This is due 
to the flake being forced to cross high potential-energy barri- 
ers. Fig. |3b), with little or no possibility to avoid them by (i) 
shifting away from the pulling direction, (ii) deforming, and 
(iii) "rotating around" {\<Pa\ < 0.2°). A similar behavior is ob- 
served also in different geometries. The spring strength tuning 
the coupling between the flake and the AFM tip plays therefore 
an important role in the model calculations. We checked that up 
to spring constants of an intermediate value {K = 0.5 eV/A^) 
our model still performs a stick-slip motion similar to the one 
found in experiment [4] (where the cantilever harmonic con- 
stants values was estimated K ^ 0.36 eV /A^), especially in the 
fully commensurate (p = 0° stacking. 



3.3 Flake-size effects 

Figure [To] compares the frictional behavior of flakes of differ- 
ent size, showing that friction tends to decrease with increasing 
flake size. In detail, the static friction force of the 3 flakes is 
Ffric = 2.99 nN, 1.38 nN, and 1.24 nN for the 24, 54, and 96- 
atom flakes respectively. This decrease is not surprising, since 
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24 atom flake 
54 atom flake 
96 atom flake 




tip position x^.^ [nm] 
Fig. 10. (Color online) The parallel component of springs lateral force, 



as a function of the support position x,ip for three flake sizes: 24- 



atom flake (solid), 54-atom flake atoms (dashed), and 96-atom flake 
atoms (dot-dashed). The three simulations involve a total applied load 
Fjv = 100 nN, 6 = 0°, = 30°, and soft springs constants K = 
0.1 eV/A^. 
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Fig. 11. (Color online) The instantaneous rotation angle (p^ as a func- 
tion of the support position jctip for three values of the support stacking 
angle: (a) if) = 0°, (b) = 15°, and (c) (p = 30°. Three flake sizes 
are considered: A'" = 24 (black solid line). A'" = 54 (red dashed line) 
and A'" = 96 (blue dot-dashed line). Simulations are carried out with 
total applied load of 100 nN, pulling angle 6 = 0° and soft springs 
a: = 0.1 eV/A-. 



the more reactive atoms at the flake boundary tend to bend 
down toward the substrate. As a result, friction is dominated 
by these boundary atoms, which amount to 75% of the 24-atom 
flake but only 44% of the 96-atom flake. Moreover, the 96-atom 
flake advances continuously, and shows no stick-slip, at vari- 
ance with the 24-atom and 54-atom flakes. This difference is 
due to a reduced rotational freedom of the 96-atom flake flake, 
due to coupling to the support acting at a larger distance from 
the flake center. The flake rotational freedom, i.e. the angular 
range of (pA explored around (f>, does represent a key issue in the 
friction physics of carbon flakes sliding over a graphite surface, 
as pointed out by Filippov et al. |[T9| . 

In the framework of our model, with each atom tied to the 
moving support by an individual spring, the flake can both shift 
normally to the pulling direction and rotate around its center of 
mass: these degrees of freedom (and, more weakly, the possi- 
bility of the flake to distort) affect friction in two very differ- 
ent manners depending on the contact being commensurate or 
incommensurate. When the flake is pulled at a commensurate 
stacking, e.g. (f> = 0°, it encounters high potential-energy barri- 
ers thus high friction: the combined possibility of rotations and 
lateral shifts allows the flake dribble the high barriers through 
local changes of trajectory. The first effect of flake shifts and 
rotations is then to reduce the friction of highly commensurate 
contacts. In contrast, when the flake slides with an incommen- 
surate stacking, e.g. = 15°, it does not encounter high-energy 
barriers nor efficiently binding configurations, thus producing 
a low-friction motion. However rotations and shifts allow the 
flake to locate deeper energy wells (both moving apart from the 
pulling direction and rotating so as to explore different stacking 
configurations), where the flake can stick, eventually providing 
sizeable friction. This second effect is therefore to raise the fric- 
tion for incommensurate contacts and eventually destroy super- 
lubricity, as was observed and discussed by Filippov et al. ||T9l . 



In our simple model the single parameter K tunes the flake rota- 
tional freedom and that of shifting perpendicular to the pulling 
direction. However, while its effect on the translational free- 
dom is independent of size, the rotational freedom does depend 
on the larger torque that springs of the same stiffness exert on 
flake atoms more remote from the flake center, as is also to be 
expected for a flake sticking to a not too sharp AFM tip. 

This point shows clearly in Fig.[TTl which displays the evo- 
lution of the actual rotation angle (pA along the scanline, for 
three different support stacking angles and three different flake 
sizes: A^'' = 24, 54, and 96. Rotational fluctuations decrease 
as the flake size increases. Indeed significant systematic de- 
viations of (pA from (p are apparent in many cases, especially 
A^'^ = 24. In particular, the small 24-atom flake for cp - 30° 
rotates all the way to (pA ^ 10°, thus displaying angular os- 
cillations in excess of 15°, for an average angle {(Pa) - 16°, 
very different from (p. When plotting the dependence of fric- 
tion force Ffiic on the stacking angle, it will make more sense 
to use, instead of the initial stacking angle (p, the average flake 
rotation angle {<Pa}, although even this indicator does not ac- 
count properly for rotational fluctuations. This large rotational 
freedom, for hard springs, is almost completely frozen: in that 
case the largest rotational fluctuation we observe is as little as 
^3°forA^'^ =24 and = 15°. 

As for lateral shifts, for soft springs K - O.l eV/A^ we 
observe shifts perpendicular to the pulling line of the order of 
1 A, depicted in Fig. [12] which displays the actual path fol- 
lowed by the center of mass of a 24-atom flake pulled by the 
support along a 6* = 0° scanline. Note that, similarly to rota- 
tional fluctuations, the incommensurate stacking (p - 30° yields 
larger lateral shifts than the commensurate stacking <p - 0°. 
With hard springs the possibility of the flake to perform lateral 
shifts is strongly reduced, so that the actual trajectory of the 
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' 24 atom flake 
54 atom flake 
■ — « 96 atom flake 



1 1.2 
X [nm] 

Fig. 12. (Color online) Subsequent points marking the trajectories of 
the center of mass of a 24-atom flake in the x-y plane for (f> = 30° 
(incommensurate stacking, blue circles) and for = 0° (commen- 
surate stacking, black diamonds) at the end of each relaxation cycle. 
The dashed lines represent the support scanlines {Q = 0°); large cir- 
cles represent the substrate atomic positions. The simulations are the 
soft-spring ones of Fig.[9l 



flake center of mass remains very close to the support scanline. 
For example, springs of K = 2.5 eV/A- yield perpendicular 
shifts < 0.2 A along the same trajectory. The rotational free- 
dom plus the lateral shifts of the flake can lead to effectively 
commensurate contacts even for an incommensurate stacking, 
thus explaining the deep energy valleys of the soft-spring pat- 
tern of Fig.|9}3, eventually responsible for the stick-slip motion 
demonstrated by the lateral force patterns of Fig. |9^. For hard 
springs constants, locking into deep energy minima does not 
occur, but at the same time the flake is driven into highly repul- 
sive geometries which it cannot dribble. This leads to higher 
force peaks, and eventually to a larger static friction. 



3.4 Angular dependence of friction 

As a reference benchmark we consider a nearly rigid flake model, 
where atoms are allowed to relax only in the z direction, corre- 
sponding to the — > oo Umit of the model studied until here, 
and comparing directly to the model used by Verhoeven et al. 
ifm . The complete suppression of angular fluctuations should 
produce an extremely sharp angular dependency of the friction 
force. 

Figure [13] shows the computed static friction force as a 
function of the fixed stacking angle (p for different flakes. We 
note that friction decreases with the angle, showing a maxi- 
mum peak centered at (p - 0°, similar to the outcome of pre- 
vious model calculation ifTTll . As shown by experiments |4|5), 
friction is maximum at an highly commensurate contact (0 = 
0°) and decreases rapidly as the flake rotates to incommensu- 
rate stackings. The friction peak is sharper for wider flakes. The 
sharpest peak for the 96-atom flake is similar to the one exhib- 
ited by the rigid model of Ref . ifTTI . Ff^c decreases by nearly 
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stacking angle (|) [degree] 



30 



Fig. 13. (Color online) Friction force Ffnc as a function of the fixed 
stacking angle (p for three different flake sizes: 24 atoms (black solid 
line), 54 atoms (red dashed line) and 96 atoms (blue dot-dashed line). 
The simulations are carried out with pulling angle 9 = 0° and total 
applied load of 100 nN. Flake atoms are allowed to relax only along z 
direction. 



one order of magnitude from the high-friction = 0° com- 
mensurate angle to the low-friction cp ^ 30° incommensurate 
one. This drop is smaller than was found by experiment lU, 
where it exceeded significantly one order of magnitude. Also 
the absolute values of friction are systematically larger than ex- 
periment. Experiment shows friction peak values near 0.2 nN, 
while the present model yields a peak value of order 10 nN, 50 
times larger This difference is even larger in the "superlubric" 
region near <p = 30°. These differences are to be ascribed to the 
larger load, the adopted short-ranged TB parameterization, and 
the neglect of thermal fluctuations, as discussed below. 

The rigid-flake models, studied here and in previous work 
[[U do not look not especially realistic, since in practice a car- 
bon flake does deform, shift and rotate while interacting with 
the graphite substrate and the AFM tip. Figure [14] reports the 
dependence of friction force Ffric on the average rotation angle 
{(Pa), for a flake whose atoms are allowed to relax in all direc- 
tions, for soft spring constants K - Q.l eV/A^. In all calcula- 
tions except those of 96-atom flake we use the same angles <p = 
0°, 5°, 10°, 15°, 20°, 25° and 30°, but the possibility of flake 
rotations allowed by the soft springs produces significantly dif- 
ferent effective average angles {<Pa}, especially for the 24-atom 
flake. Not surprisingly, with its vast rotational freedom, the 24- 
atom flake displays an almost ^-independent, constant friction. 
For such a small flake with soft springs, rotations and shifts are 
so effective to hinder the possibility to observe any reliable (f>- 
dependency of Ffnc- The 54-atom flake and, more clearly, the 
96-atom flake show average angles nearer to the support val- 
ues, with smaller-amplitude fluctuations, and therefore display 
a friction curve behavior with a peak at {<Pa} - 0°, fairly sim- 
ilar to the one obtained in the semi-rigid case and observed in 
experiment, and with smaller friction at incommensurate an- 
gles. These results suggest that when the FFM tip happens to 
bind to a graphene flake constituted by substantially less than 
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Fig. 14. (Color online) The friction force Ffn^ as a function of the aver- 
age rotation angle (0a> for three different flake sizes: 24 atoms (black 
solid line), 54 atoms (red dashed line) and 96 atoms (blue points). Sim- 
ulations involve a pulling angle 9 = 0°, total applied load of 100 nN 
and soft springs constants = 0. 1 eV/A^. 
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Fig. 16. (Color online) Nonlinear dependency of (a) the friction force 
and (b) the friction coefficient [i = on the total applied load Ffj , 
for a commensurate contact (ip = 0°, black solid line) and an incom- 
mensurate contact (0 = 15°, red dashed line) of the 24-atom flake 
flake. Simulations are carried out for pulling angle = 0° and rigid 
springs, K = 2.5 eV/A^. 



ing that for realistically large flakes in excess of one hundred 
atoms, the precise value of the spring stiffness should become 
irrelevant, as long as it remains in the < 1 eV/A^ region. 
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Fig. 15. (Color online) Comparison of the friction force dependence 
on the average rotation angle (0^) for a 54-atom flake for two values 
of springs strength: soft A" = 0.1 eV/A^ (black solid line) and hard 
K = 2.5 eV/A^ (red dashed line), plus the semi-rigid case (blue dot- 
dashed line). Simulations involve a pulling angle 9 = 0°, total applied 
load of 100 nN. 

approximately 10^ atoms, no clear angular dependency and no 
superlubric regimes are observed. 

Figure [15] summarizes the effect of increasing the rigidity 
of the springs on the friction dependency on {(Pa)- the fric- 
tion peak at a commensurate arrangement becomes sharper and 
sharper as the spring rigidity increases. At variance with the 
radical changes in {cf>A} dependency of the 24-atom flake as a 
function of the spring constant, for the 54-atom flake the shift- 
rotational effects become comparably less important, suggest- 



3.5 Load dependency 

We come now to study the dependence of the friction force Ffac 
on the applied load F^, exploring a range 20 to 100 nN, match- 
ing typical experiment values II4I29I . Figure [T6l shows the de- 
pendence of the friction force Ffi-ic and coefficient /i s Ffnc/F^ 
on the applied load. Hard springs are selected to reduce the 
flake shift-rotational effects, in order to focus on the load de- 
pendence of friction and simpler comparison with earlier re- 
sults. Friction increases with load, but significant deviations 
from the linear Coulomb law are observed, especially for com- 
mensurate stacking = 0°. Observe that experiment found an 
even weaker dependency of the friction force on load [:4!]. Al- 
though the data do not point clearly in the direction of a power- 
law behavior Ffi-ic <^ F"-^, it is clear that if any such law was to be 
estimated, it would have a < 1. This is at variance with previ- 
ous findings for a sharp undeformable tip-surface contact ITSl . 
and with recent studies of the sliding of hydrogen-passivated 
carbon 1301 . The resulting static friction coefficient approaches 
the standard macroscopic value 1311 of graphite-graphite con- 
tact, yu = 0.1, while much smaller values are found in the single- 
crystal FFM experiments addressed by the present model. Re- 
gardless of the applied load, the flake-substrate distance being 
smaller in the model than in real life produces larger absolute 
values of friction, and this overestimation is particularly severe 
at small load. 



3.6 Scanline dependency 

We investigate the dependence of the friction force versus stack- 
ing angle on the actual scanline followed by the support. Chang- 
ing scanline determines a different potential profile seen by the 
flake, thus modifying the frictional behavior IITSl . Figure [TTlre- 
ports the friction angular dependency for three equally spaced 
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Fig. 17. (Color online) (a) Friction force Ffdc as a function of the 
average stacking angle {0/i> for three different scanlines drawn in 
panel (b), defined by the three following initial stackings of the sup- 
port over the substrate: AB (black solid), AB with a transverse shift 
h = <igraph/4 (red dashed) and AB with transverse shift h = dgraph/2 
(blue dot-dashed). The simulations involve a 24-atom flake, pulling 
angle 6 = 0°, applied load Ffj = 100 nN and hard springs of constants 
K = 2.5 eV/A-. 



scanlines. We carry out simulations for hard springs K -2.5 eV/A 
where scanline effects are the most visible, because of hin- 
dered lateral shifts. For the h - Q and h - (igraph/2 scanlines 
we find a similar friction for all values of the stacking angle 
{<Pa), while the h — iigraph/4 scanline shows systematically 
lower friction, especially for small 0. The reason is that for 
4>A — Q along this special line each flake atom never hits any 
substrate atom directly on top, therefore effectively finding a 
significantly lower corrugation. In contrast along the two other 
scanlines, for = 0° one half of flake atoms encounters peri- 
odically a substrate atom right below its trajectory, thus finding 
a high corrugation. Softer springs produce a much weaker de- 
pendence on the scanline: the flake takes advantage of its free- 
dom to displace laterally, thus following low-corrugation lines 
(such as h - c/graph/4) even when the support pulls it along 
some nearby parallel line. 

The scanlines of Fig. [17] involves - 0°, i.e. a pulling 
along the x direction, where the flake encounters periodic rep- 
etitions of the substrate potential. A different pulling angle af- 
fects directly this periodicity of the problem, in general leading 
to a nonperiodic profile. Figure [18] displays the friction force 
as a function of the pulling angle 9. Data are averaged over 
three different scanlines, defined by initial stacking with h - 0, 
h - (igraph/4 and h - (igiaph/2. Like in previous calculations 
lim . we find a minimum friction for pulling angle 9 - 0°, fol- 
lowed by a fast growth in friction (until 9 - 10°). We attribute 
the observed differences between our results and those by Ver- 
hoeven et al. to the diff'erent interaction models. 



Fig. 18. Friction force Fine as a function of the pulling angle 6. Static 
friction data are obtained by averaging three different scanlines, de- 
fined by initial stacking AB, AB shifted by rfgiaph/4 and AB shifted by 
'^eraph/2 perpendicular to the pulling direction. Simulations involve a 
24-atom flake with load of 100 nN, support stacking angle (p = Q° and 
springs constants K = 2.5 eV/A^. 



4 Discussion and conclusion 

We find fair qualitative agreement between the results obtain- 
ing by our TB atomistic model and the existing experimental 
data, with a few significant difi'erences. Firstly, our calculations 
recover the stick-slip behavior of the lateral forces, character- 
istic of FFM sliding experiments. In particular we find the cor- 
rect qualitative dependence of stick-slip on the springs stiffness 
characterizing the cantilever-tip-flake coupling: soft springs al- 
low for clean stick-slip behavior while hard springs inhibit it. 
Our calculations also reproduce correctly the friction pattern as 
a function of the average stacking angle {<Pa} especially as long 
as the rotational degree of freedom (pA is quenched. We also 
find that for larger flakes, the fluctuation in <pA are suppressed 
automatically anyway, due to the larger torque exerted by the 
springs connecting the flake to the tip. Accordingly, for flakes 
of sufficiently large size incommensurability produces signifi- 
cantly less friction, although the friction drop is smaller than in 
experiment. In the quantitative comparison between the exper- 
imental results and our model, we find static friction force Ffric 
and coefficient fi systematically at least one order of magnitude 
larger than experiment, this difference being especially signif- 
icant in the incommensurate configurations where no proper 
superlubric regime is observed. These and other quantitative 
discrepancies are to be attributed to: (i) The reduced interlayer 
equilibrium distance, related to the small cut-off' distance of 
the present TB parameterization, which is responsible for the 
increased energy coiTugation experienced by our model flake 
with respect to real graphene on graphite, (ii) The extra reac- 
tivity of the isolated model graphene flake with respect to a 
real one, which is bond to the AFM tip and thus somewhat 
passivated; accordingly, especially the atoms at the flake bor- 
der, show a greater tendency to react with substrate atoms, 
thus increasing friction, (iii) The neglect of thermally-activated 
slips through energy barriers 0271281 : this neglect generates 
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an overestimation of friction especially where these barriers 
are lower, i.e. at incommensurate stackings. Indeed the cur- 
rent understanding |32| of the observed [:4| superlubric slid- 
ing involves thermolubricity associated to a high attempt rate 
for overcoming the corrugation barriers due to the microscopic 
mass of the vibrating tip apex. If the experiment of Ref. ^ 
could be repeated at the much lower temperature of a few de- 
gree Kelvin, the observed friction values and dependency on 
the angle would probably look much more similar to the one 
obtained in the present model. 

Calculations carried out with comparably soft springs and 
small flakes {Na < 54) show that the flake shift-rotational free- 
dom increases friction for incommensurate stackings (by al- 
lowing the flake to explore deeper-bound minima) and decreases 
it for commensurate ones (by allowing the flake to dribble the 
top barriers): the result is a substantial flattening of the de- 
pendency of the friction static force Ff^c on the stacking an- 
gle (p. Harder springs (e.g. K = 2.5 eV/A^) would suppress 
the flake freedom to rotate and shift laterally but are incom- 
patible with the clear stick-slip behavior observed experimen- 
tally. These considerations confirm that the size of the flakes 
showing superlubric sliding in actual FFM experiments is large 
A^fi > 96. Many discrepancies with experiment would proba- 
bly be disposed of, if a longer ranged interatomic interaction 
was employed, for example a TB model based on a longer 
cutoff lf33l . This way, a much weaker flake-surface interaction 
would effectively correspond to comparably stronger tip-flake 
interaction, thus a significant hindering of rotations and trans- 
lations even with a realistically weak tip-flake coupling K < 
0.5 eV/A^. If accurate long-range C-C interactions up to dis- 
tances of the order of 1 nm were present in the force field, one 
could even include substrate relaxation to study an even more 
realistic model. Such an improved model would however im- 
ply significantly larger computational workload, especially if 
thermal excitations were also included. Further research should 
also investigate the effect of the presence of structural defects 
in the flake or in the substrate, as proposed by Guo et al. Il34l . 
and the effect of flake shape. 
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